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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE 23^ 

A SURVEY OF METHODS FOR DETERMINING STABILITY 
PARAMETERS OF AN AIRPLANE FROM DYNAMIC 
FLIGHT MEASUREMENTS 
By Harry Greenberg 


SUMMARY 


Various methods of reducing to stability parameter form the 
response to sinusoidal and transient disturbances are discussed, using 
the simplified longitudinal motion of an idealized airplane as an illus- 
trative example. It is shown that there are basic limitations in the 
determination of some of the stability derivatives as compared with the 
trans f er~func t i on coefficients, which are certain combinations of sta- 
bility derivatives directly related to the airplane response. Hence, 
most of the report is concerned with methods of determining transfer- 
function coefficients rather than stability derivatives. 

It is shown how the method of least squares can be applied to give 
the desired parameters and also the ratio of their error to that of the 
basic data. The determination of these parameters and their correspond- 
ing error ratios is a nonlinear problem which it is shown can be solved 
by linearization using a first approximation to the unknown parameters. 

A number of methods of obtaining a good first approximation, which also 
involve a least squares procedure, are explained and illustrated in the 
numerical examples. 

Although the examples are confined to a simplified case of longi- 
tudinal motion, the methods presented are applicable in general to other 
more complicated types of motion. 


INTRODUCTION 


The scarcity of reliable data on the stability characteristics of 
aircraft at transonic and supersonic speeds and the difficulties of 
obta inin g this information from wind-tunnel tests (particularly with 
regard to dy nami c parameters, such as the rotary damping derivatives), 
have accelerated interest in methods of obtaining such data from flight 
tests. Also the extensive use of automatic stabilization and . control 
equipment and the uncertain and generally poorer dynamic-stability 
characteristics arising from the use of unconventional configurations 
necessitate comprehensive and refined measurements of the stability 
derivatives . 
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In the past, limited information has been derived from measure- 
ments of airplane characteristics in steady-straight and steady— turning 
flight. Recently, frequency— response measurements to sinusoidal control 
deflections have been employed to evaluate additional dynamic— stability 
parameters. The brief testing period available during flights of 
missiles and high-speed research airplanes has impelled consideration of 
transient— response flight— testing methods and associated instrumenta- 
tion. Proper testing and analysis techniques should permit evaluation 
of all the stability derivatives, including the limited number which can 
be evaluated from steady— flight data. 

The oscillation technique was discussed by Laitone, Cornell Aero- 
nautical Laboratory, in an unpublished report on dynamic flight measure- 
ments in which the response of the airplane to sinusoidal elevator 
deflections was measured. It was shown that the lift-curve slope, 
elevator effectiveness, and damping and stiffness parameters associated 
with the short-period oscillation could be determined. A more detailed 
presentation of this method and an application to flight data are given 
in reference 1; similar theoretical and flight-test investigations for 
the lateral motions have also been reported. The method of analysis of 
reference 1 is only applicable to simple dynamical systems — that is, 
systems which mathematically are similar to one with a single degree of 
freedom, and is incapable of reducing the stability parameters to the 
basic stability derivative form. 

The obvious advantage, from the standpoint of test simplicity and 
time, of using the response to a step elevator input instead of the 
frequency-response tests was soon realized, and the work reported in 
reference 2 shows how the step-response data can be converted into the 
frequency— response form. Subsequent work (reference 3) extended this 
method to the response to a pulse elevator input. Another method of 
analysis of response to an arbitrary elevator motion was suggested by 
Loring and Jonah of Chance— V ought Aircraft Company. This method does 
not require transference from the "time domain" to the "frequency 
domain" as is the case with references 1 through 3j and- is referred to 
later in this report as the "derivative method." In this derivative 
method there appears for the first time the application of the method 
of least squares to obtain the most reasonable values of airplane 
parameters from redundant measurements . 

Examination of the available literature indicates a lack of infor- 
mation concerning the general methods of analysis applicable to more 
complicated systems , for example, systems with more degrees of freedom 
or with higher-order derivatives. The purpose of the present report is 
to establish more general and rigorous methods for determining aero- 
dynamic parameters from dynamic flight measurements. The following 
principal and basic problems are studied: the relation between the 

number and type of applied forcing functions and measured responses and 
the corresponding number and type of determinable aerodynamic parameters 
various methods of converting flight data to a form suitable for 
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determination of aerodynamic parameters; and the correct application of 
the method of least squares to compute the aerodynamic parameters. 

Although the methods presented sure applicable to more complicated 
systems, the examples in the present report are confined to the simpli- 
fied longitudinal motions of an idealized airplane in order to facili- 
tate computations. 


NOTATION 


General 


a angle of attack, radians 

angle of attack of tail, radians 
A,B in- and out-of— phase components of oscillation 


D 

8 

€ 

E 

E r 

El 

g 

0 

% 

X 

L 

La 

LS 


differential operator 

elevator deflection, radians 
dovnwash angle, radians 
residual error in an equation 
residual error in a real equation 
residual error in an imaginary equation 
acceleration due to gravity 
angle of pitch, radians 

pitching moment of inertia, slug— feet squared 

roots of characteristic stability equation 

lift force, pounds 

dL 

da 

dl 

dS 
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k 


Lfc 


Iry 

m 

M 


M, 


a 


M S 


Mg 

My 

n 
< 1 

R,cp 

S 

t 

V 

Vo 

¥ 

w q 

x 


CD 


tail lift, pounds 
9l 


mass of airplane, slugs 
pitching moment, foot-pounds 


9m 

9a 

Jim 

9tDa 

9m 

35 

9m 

9q 

9M 

9v 

relative density parameter 



2 x mass 

density X wing area X wing chord 


normal acceleration, g units 

angular pitching velocity, radians per second 

amplitude and phase of complex number 

sum of weighted squares of residuals 

time, seconds 

velocity of airplane, feet per second 
trim velocity of airplane, feet per second 
weight of airplane, pounds 

weighting factor indicating accuracy of q measurement, etc. 

longitudinal distance between center of gravity and neutral 
point of airplane, feet 

angular frequency, radians per second 
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Subscripts 


c calculated 

m measured 


Transfer-Function Coefficients 


b 


damping parameter 


f “Mq-MDa \ 

V Iy mVo ) 





-M. 

m Q 

MDa 

Iy mV 0 Iy 

Ms Ls Mrt 
!y + mV 0 Iy 


C M 5 L q Mq L S 

°9 Iy mV Q Iy mV Q 

°On °Oq ^ 

(Cl q ~ Co a ) ^ 

°2n ~ c ia if 


k 


stiffness parameter 



h, 

mV 0 I y J 


MATHEMATICAL AND AERODYNAMIC PRELIMINARIES 


In order to illustrate the types of dynamic parameters involved 
and the conditions under which, they may be measured, the case of longi- 
tudinal motion of an airplane will be considered. In setting up the 
equations of motion, the following assumptions are made: 
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Basic Assumptions 


1. Linear equations with constant coefficients 

2. Measurement of forcing function not subject to error 

3. Response measurements subject only to random errors 
Constant airspeed and level flight 

5* Aerodynamic lift equals al^ + SLg 

6. Aerodynamic moment equals oMq, + Da Moa, + q Mq + 5Mg 

7* Rigid airplane 

8. No unsteady lift effects except for downvash lag -which 

introduces the derivative Mpa. Only the first three of the 
above assumptions are necessary in order to apply the methods 
of this report; the others are made to simplify the numerical 
examples given later and because they do approximately 
describe the airplane motions in the maneuvers which are con- 
sidered in this report. 


Equations of Motion and Statement of "Inverse" 
Problem of Airplane Dynamics 


Based on the above assumptions, the longitudinal equations of 
motion may be written 


— ala - Da mVo + <1 mV 0 = L 


( 1 ) 


— aMct, — Da Mp<x — qMq + IyDq=M 


( 2 ) 


The forcing functions, lift L and moment M, may be applied aero- 
dynamically by deflecting an elevator, as discussed in a later section 
of the report, or nonaerodynamically as, for example, by releasing a 
bomb or firing a gun. 

Most of the applications of the mathematical theory of dynamic 
behavior of airplanes have consisted of computations of airplane 
response characteristics when the stability derivatives such as 1$,, 

Mo,, and Mq are assumed known. The problem considered in the present 
report is the inverse; namely, given the airplane response in a, q, 
or n to a disturbance, to evaluate the stability derivatives of the 
airplane . 
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It will be proved later in this report (see section "Discussion") 
that the determination of moment derivatives from flight data is subject 
to certain basic theoretical limitations. There are certain combina- 
tions of moment derivatives, however, which determine the behavior of 
the airplane and which can be computed from the flight data. These 
combinations of derivatives are called "transfer coefficients" and are 
defined below. Most of this report is devoted to the determination of 
these transfer coefficients. 


Transfer Functions for Control Deflection Input 


For this case, in which L=SL§ 
tion of equations (l) and (2) is 


and M=5M$, the operational solu- 


D5 ^5 ^8 Mq 

+ -y + ^-Q-y 

|2+D f_ + ^ - 

\ Iy Iy mV 0 / mVo Iy Iy 


(3) 


D 2 +bD+k 


f Mg _ _Ls_ Mpq N % MS bq 

q Vly nV Q I 7 ) D ~ Iy mV 0 + I y mV Q 

8 D 2 +bD+k 


(5) 


Ci q D+C 0q 

D 2 +bD+k 


( 6 ) 


Since the normal acceleration is often measured, it will be useful to 
give the solution 

& _ Vo Dec _ Vo iqP a+ ( C iq~Doc{,)D + Doq 

5 g 5 g D 2 +bD+k 


C2 n D 2 +Ci n D+Co n 


D 2 +bD+k 


( 7 ) 
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Quantities such as Co a , Ci a , C 0 q, Ciq, b, and k will be called 
"transfer coefficients," because they determine the transfer function 
between a response and a disturbance. 

It may be noted that there is a direct relation between the transfer 
function and frequency response. The sinusoidal response of the airplane 
to a sinusoidal disturbance of frequency us is obtained by substituting 
ico for D in the corresponding transfer function. The magnitude and 
phase ang le of the complex number so obtained give the amplitude ratio 
and phase (lead) between the response and disturbance, which are of the 
same frequency. 


COMPUTATION OF DYNAMIC PARAMETERS 
(TRANSFER COEFFICIENTS) FROM 
FLIGHT DATA 


A more detailed study of the analysis procedures for reducing basic 
data obtained by various flight— test methods is now of interest. As 
previously noted, a dynamic flight test made for the purpose of measur- 
ing stability and control characteristics generally consists of the 
measurement of an input disturbance (such as control deflection) and the 
resulting airplane response in one or more degrees of freedom. The dis- 
turbance and response may be transient or may be sinusoidal motions, 
with time and angular frequency, respectively, as the independent 
variable. Equations of motion are then assumed, and the associated 
transfer coefficients and (to the extent it is possible) the stability 
derivatives are evaluated by a process usually referred to as "curve 
fitting." The mini mum number of points required to fit a curve to the 
observed data is the same, of course, as the number of assumed transfer 
coefficients. A larger number of points ordinarily are obtained with 
either transient or sinusoidal tests, and the values of the stability 
parameters which will fit these redundant data with the least error 
should be determined. This is done by application of the principle of 
least squares, which also furnishes information about the accuracy of 
the parameter evaluations in terms of the accuracy of the basic data 
( error ratio) . 


Principle of Least Squares 


This principle states that the most probable values of the unknown 
parameters are those for which the sum of the squares of the errors 
(differences between calculated and measured airplane motion, for 
instance) is a minimum. If only one airplane response is measured, for 
example, the transient pitching velocity q from t=0 to T, then the 
inverse problem of airplane dynamics is solved by finding values of 
b, k, Ciq and Coq such that 



NACA TN 2340 


9 


S = 



2 

Unf-Sc) dt 


is a minimum. 

If more than one airplane response is measured, for example, the 
transient responses in 0, q’, and Dq from t=0 to T, then the values 
of the parameters h, k, Cig and Cog are to he determined such that 


S 



(0j)f-0 C ) 2 dt + Wq 



(qnr^c ) 2 dt 



is a minimum. The weight w of a measurement is a number indicating 
the accuracy of that measurement (as regards random errors). More 
specifically, the weight is the reciprocal of the mean square error 
(normal distribution in the errors of measurement assumed, see refer- 
ence 4) . 

If the frequency response (both amplitude and phase) is measured 
over a range of frequencies 03, then the parameters are to be deten- 
mined by the condition that 

S = y 1 wr (RuH^c) 2 + vcptcpur-cpc) 

is a minimum where the summation is taken over the frequencies at which 
the response is measured. 

The problems above are nonlinear in the unknowns b, k, Ci, and C 0 . 
The only practical method of solution is to linearize the problem and 
iterate from a first approximation. lor the case in which only one 
quantity is subject to error (like the first case mentioned above) the 
methods of linearization and iteration are explained in reference 4, 
pages 214 and 84, respectively. The method of linearization is dis- 
cussed briefly in the next section and in the appendix. The subsequent 
section (which constitutes the principal part of this report) deals 
with methods, most of which also involve least squares solutions, 
of obtaining a good first approximation. 

The idea of determining the parameters from a transient response by 
linearization and iteration from a first approximation is due to Shinbrot 
of Ames Laboratory and is discussed and exploited more fully in refer- 
ence 5* 
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Deter miration of the Parameters and Their 
Relative Accuracy by Linearization 


To determine the parameters b, k, Ciq, and Coq from q and 5 
measurements, it is first necessary to determine an initial approxima- 
tion (b 0 , kQ, Ciq Q , and Coq Q ) to the parameters (e.g., by one of the 

methods of the next section). Then linearization of the problem of 
deter mining the incitements in the parameters vhich make S a minimum 
gives the equations: 


where 


and 


4 

^ Sij xj = Si 1 = 1, 2, 3, 4 
i=l 


s ij 


d g S 

^Xi^Xj 


Si 


dS 

3xi 


Xi = b-h 0 x 2 = k— ko 


** * %■"%„ 


The partial derivatives are to be computed at the values b=b Q , k=ko, 

Ci =Ci and C 0( _=C 0 _ . Expressions for Sij and Si for the transient 
h 1o ’ "o 

response and frequency response are given in the appendix. 


Having obtained the increments X*, it may be necessary to repeat 
the calculation using the corrected values of the parameters as a new 
first approximation and so on. This, of course, will not be necessary 
if a good first approximation is obtained. It appears to be worth 
while, therefore, to devote some effort to obtaining a good first approx- 
imation. In the next section it. is shown how the method of least squares 
is used to find such an approximation. 

The relative accuracy or weight (reciprocal of mean square error) 
of a computed parameter is obtained as shown in reference 4 from the 
following formula: 


1 cofactor of Sii 

^I = fsTJ] 
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A more useful indication of the accuracy of a computed parameter 
for some purposes is the ratio of percentage error in the computed 
parameter to the percentage error of the flight measurement. This ratio, 
which may he called the error ratio, is (for q measurements) 


percent error in 
percent error in q 


= error ratio of = 





cofactor of Sj.i 

|Sij| 


where is the full-scale reading of the q instrument, and w q 

is the weight of the q measurements (reciprocal of mean square error), 


Methods for Obtaining a First Approximation 
to the Parameters 


Sinusoidal (frequency) response .— In this case the measurements 
consist of input and response amplitudes and phase shifts between the 
input and response in one or more degrees of freedom over a range of 
frequencies. For example, in the case of longitudinal motion with two 
degrees of freedom the transfer function between pitching velocity q 
and elevator angle 5 has been shown to be: 

I = Cl q D+G °q 

5 D 2 +bD+k 
If 

5 = 5 0 sin cut 


and 

q 

= R sin (cot+cp) = A sin cot + B cos cot 
D o 


then it can be shown that A+iB is the value of the transfer function 
when D=ico. This gives the equation: 


A(k-co 2 ) — Bbco + i[B(k-co e )+Abco] = Co q + Ci q coi 
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Equating real and imaginary parts gives the two equations: 


Ah — Baib — Coq = Ato 2 Real equation 

Bk + Aaib — Cig to = Bai 2 Imaginary equation 

If the values of A and B are known at two frequencies, then it 
is possible to set up four simultaneous equations in the four unknowns 
b, k, Coq, and Ciq as follows: 


Aik — BiCOib — Coq = Ai03i 2 

A 2 k - B^b - Coq = A2CO2 2 


Bik + Aioiib 


“ ^lCiq = BitDi 2 


B 2 k + Asfl^b 


“ O^C = BgOfe 2 


where Ai, A 2 , B 1# B 2 are values of A and B at asi and 032 . 

The solution of these equations gives values of b, k, Ciq, and Coq 
corresponding to a frequency response curve that agrees exactly with the 
measured curve at the two frequencies selected. This has been called the 
"method of selected points." 


A better first approximation can be obtained at the expense of some 
additional labor by using data obtained at more than two frequencies. 

The method of least squares is used to calculate values of the parameters 
that make ZEr 2 + ZEp 2 a minimum where Er and Ep are defined by: 


Ak — Bo3b — C 0 q — Acs 2 = Er 
B k + Amb — Ciq 03 — Ben 2 = Ep 
This leads to the following set of normal equations: 

2 ( B 2 03 2 +A a 032 )b + Z Bo3C 0q - 2 A^C i q = 0 

Z(B 2 +A 2 )k -ZAC 0 q - ZBcuCiq = Z(B 2 a32+A 2 a)2) 

ZBodb -Z Ak + NC 0q = -ZAm 2 

ZAa^b +ZBak -Za^Ciq = ZBcd 3 

where E equals the number of frequencies and the summation is taken 
over all the frequencies at which data are obtained. 
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If the phase and amplitude of the normal acceleration response to 
sinusoidal elevator motion is measured, the coefficients of the transfer 
function 


n = Ca^+CinD+Con 
8 D s +bD+k 


can he determined in a similar manner. 

To illustrate the computation of a first approximation hy the 
method of least squares, the frequency response of an idealized airplane 
to sinusoidal elevator motion was calculated and the results were used 
to compute the transfer coefficients. The frequency— response points 
used in the computation are shown in figure 1. Table I gives a complete 
schedule of calculations of the transfer coefficients k, b, Ciq, and 
C 0 q from the pitching— response data of figure 1. 

Transient response .— Since it is much faster and easier to measure 
the transient response than the frequency response of a dynamical 
system, the analysis of transient responses to determine transfer coef- 
ficients is of great interest. A common type of transient function, the 
pulse, is shown in figure 2; some special types of input functions 
(impulse, step, and ramp) are shown in figure 3. Several methods of 
analysis, of varying degrees of generality, will be explained. 


(l) Inspection of the transient .— If 8 becomes constant after 
a brief transient period, b and k can be determined from the 
damping and period of the oscillations (assuming that the system is 

less than critically damped) . If T 1/2 is the time for the free 
oscillations to damp to half amplitude and P is the period, then 


b 


1.386 

T l/2 


and 


_ 0-48 + 39.48 

' T 1/2 2 P (l) 2 


If, in addition, the steady— state value of q and the steady-state 
value of 8 are known, then their ratio q^/Soo is equal to Coq/k. 
If the values of q^ and 8« are zero (as for a pulse input) then 



If the input is a step, then 
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Otherwise there is no way of determining Cig "by inspection. 

(2) Fourier transform .— The vector ratio of the Fourier 
integral of the output to the Fourier integral of the input gives the 
frequency response. That is, if 


/ 0 q(t)e iajt at 
/” S(t)e _i£l>t dt 


= RE ic P 


then the response to & = sin cut is q = R sin (cut-Kp). The applicar- 
tion of the Fourier integral to the conversion of transient to fre- 
quency response is discussed in references 2 and 3* 

If, as in the case of a step or ramp input, the slope of the 
output or input becomes essentially constant (and equal to Dq^) when 
t > T, the Fourier integral may be evaluated by the alternative sub- 
stitutions 

qe-ixut dt - i / T ericot dq - Dq T er^T 
o o 

or 

rp - «— iofE tv, io>T 

q e — icut ^ - 


where q^ is the value of q at t = T. 

The transfer coefficients can, of course, be obtained from the 
frequency response by the methods of the previous section. An illus- 
tration of the application of the Fourier transform to compute fre- 
quency response from the step response of figure 4 is given in 
table II. The assumed transfer function is slightly different from 
that used in table Ij hence a slightly different frequency response is 
obtained. The frequency response obtained as an intermediate step in 
the computation of transfer coefficients by this procedure is often 
itself of interest in problems of automatic stabilization. 

The so-called "incomplete Fourier transforms" qe - *^ dt and 
/q 5e — ' dt can also be used to compute the parameters from a pulse 


/ 


qe 


r lcot dt + 




-L/'irjic 


io> 


0)2 
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which has not reached a steady state and for nonzero initial condi- 
tions. Assuming the transfer function 


„ Ci D+C 0 
£ = 3 _j 

5 D 2 +bD+k 


then it can "be shown that 


[ ( ico) 2 +bia>fk ] / qe icu ^ dt + e (Dq^4hqip+iajq T ) — Dq Q — hq 0 — icuq c 
o 

= ^C lq ia>fCoq) y T Se“ iait dt + C lq 5Te“ ia)t - C lq 8 0 


from which a real and an imaginary equation can he set up. These can 
he used to set up four simultaneous equations using computations at 
two frequencies, to get approximate values for h, k, C lq , and Coq. 

If more than two frequencies are involved, then the method of least 
squares can he used as in the case of frequency— response measurements . 

(3) Derivative method .— Another method for computing the transfer 
coefficients from the transient response consists in using the 
measured values of a sufficient number of higher derivatives of input 
and response in the assumed transfer function. For example, if the 
assumed form of the transfer function is 


q _ D0 


S 5 


D 2 +bD+k 


it is possible to determine the parameters from four sets of simulta- 
neous measurements of 0, D0, D 2 0, 5, and / Q 5 dT by fitting the 
equation 


D 2 0 


+ bD0 + k0 - 


Ciq 5 — Coq 



SdT 


0 


If the subscript i denotes the value of a quantity at t=t x etc., 
then four simultaneous equations for the parameters b, k, Ci q , and 
can be set up as follows: 



H 

CJ\ 


(D 2 0)i + Td(DS)x + &9i — GiqBi 
(D20) 2 + b(D0) 2 + k0 2 - C aq 6 2 
(D 2 0)g + t(D0)g + k0g ~ ClqSg 
(D 2 0) 4 + b(D0) 4 + k0 4 - C lq 5 4 



Tlie parameters computed from the above equations can, of course, he used to compute a 0 
response for a 6 input which passes through the four given values of 5 and the integral 

of which corresponds to the four given values of /£ 6dT. However, the transient 0 curve 

so obtained will not necessarily pass through the four given values of 0, but the above 

relation between D 2 0, D0, 0, 8, and /£ 8dT will be exactly satisfied at the four points. 

If it is desired to use more data (taken at more than four time instants) then the method 
of least squares may be used to compute values of b, k, Ciq, and Coq such that £ Ei 2 is a 
minimum 

where 


Ei 


(D e 0)i + b(D0)i + k(0)i - CxqSi - Cgq 



5 dr 


The normal equations are: 


Oil£S HI VOW 



1 <M) l 2l> + i (De) i( e) i k - S(DS) 1 (B) i C lq _E(D0 ) 1 r H 8dT c 0 --E (DSJitD'eJj 

\ <DS) i (9) 1 b + l <- s h S * - c, - S(a) ± f H Mr 0 o (B« 8 ) 

1 H iJ 0 q L 1 

i (D0)i(5)ib + i (0)i(5) i k - Z (&i) 2 C lq - E (&)* y' 1 5dTC 0q = - Z(6) ± ( D 2 0) i 

? (D0> ^ !4Tb 1 (9, i{ 1 SdTk-EfS^ 1 SdTC^- C ° q =-?X tlMT 

An example of the application of the derivative method to the determination of the parameters 

??T r9 T nS9 to a step input in 6 is shown in figure 5, La L 

a 6 ^ or a step response the differential equation relating 0 nnri 6 is 

D 2 0 + b D0 + k0 — C ^ - C D t = 0 or E ’ 

1 q 

and the normal equations resulting from the minimizing of Z E 2 are 

Z (Dfl)^ b + E _e (DQ)^! -E (DO)^^ = -Z (D0)., (D 2 0) • 

i i i qi qi 11 

z (D0) i (0) + s(0) ± a k -Z(0)iG - Z0itiC o = -Z (D 2 ©) . 0 - 

1 X XuiQ.^J-1 

\ (D0) i b + \ < 0 >i k - * C i q - 2 tiC Qq o — E (D 2 0)i 

Z (D0).tb + £ 0^. k - £ t. C lq - E t t 2 C Qq = -E (D 2 0) i t 

where ET is the number of time instants. 


H 

-J 


HA.CA TIT 23^0 
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( 4 ) Prony* s method. — Frony* s method (reference 6) for fitting 
a sum of exponentials to a number of equally spaced ordinates may be 
used to obtain the transfer coefficients from a step or impulse 
response. In the case of a step it can be assumed that q can be 
represented as a constant plus the sum of V exponentials 


q = q tt + Aie^ 1 ^ + A 2 e^ 2t + . . . A v e ?tvt 

and the values of q(t) at n equidistant instants A seconds 
apart will satisfy a difference equation of the form 

q (m+v)A + &vq (m+V— l)A + * * * ai V* = 0 

where q^, q oA , . . . , q^ indicate the value of q— q^ at 
t=A, 2A, . . . , mA, respectively. The following set of 
n— V equations may be written: 


^(V+l)A + ^V^VA ■*■••• a i^-^ “ ® 
q (v+2)A + a v q (wl)A + * ’ • a i^aA = 0 


Vi + a V q ( n -l)A + ’ * * a i q (nt-v )A " 0 


If the value of q B is known, (i.e., if a practically steady 
state has been reached) then the above equations can be solved 
for the coefficients aj by the method of least squares (if n>2v). 
If the value of q« is unknown (because a steady state has not 
been reached at the end of the record), then the difference 
equation may be written, calling 

q as q m at t = nA 

V+V + ®vV+V-i + • • • a x<lia = ( a i + - * • a y + b)<leo = "^v+i 

This leads to a similar set of difference equations with the 
additional unknown a v+ r These also can be solved by the method 
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of least squares (if n>2v) for the coefficients a^, from which., 
of course, q w is determined. The least squares solution for 

n-V 

the coefficients aj_ is determined by minimizing 2 E^ 2 where 


= ^-m+v + ®v ^m+v — 1 + • • • a i1m + a V+i 


The values of are then used to form an algebraic equation 

x v + a v x v_1 + . . . a x = 0 

the roots of which are related to X (the characteristic roots) 
by 


X = 


_ log x 


Knowing the V values of X the set of n equations 
q A = A i e Xl ^ + A 2 e^^ + . . . A y e^ v ^ 


q sA = A x e 2X i A + A 2 e 2XaA + ... Ay e **V A 


^nA 


= A, e 


nXj A 


+ A-, e 


nXgA 


A v e 


n^A 


can be set up, and the method of least squares is used a second 
time to determine the unknown coefficients A^. The transfer 
function is now completely determined. To express it in terms 
of the differential operator D, each term e^ is replaced 
by D/D—X, its Laplace transform multiplied by D, and the 
resulting fractions are added, that is, 

A,D A_D AyD 

q = q + — + -S_ + . . . -k— - 

D-X 1 D-X 2 L-Xy 

Gy D V + Gy_^ D v ^ + • • • Cq 

D V + B v _ x D^ 1 + . . . B 0 
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where the coefficients are functions of the X roots, whereas 
the coefficients are functions of the coefficients Aj and of 

the X roots. 

A numerical example illustrating the application of Prony’s 
method to the response to a step input is shown in tables 17(a) 
and 17(h). This e xamp le is based on the computed response shown 
in figure 4 . The form of transfer function for the pitch response 
to elevator deflection is the same as that assumed previously, 
that is. 


1 = D + °°q. 
8 D 2 + bD + k 


The values of q are tabulated at equal intervals of time A of 
0.1 second, from t=0 to t=1.0 second, and the solution is 
carried out on the assumption that the steady— state value of q 
is unknown. 

The values of b and k are determined in table 17(a); the 
values of C lq and C 0( ^ are determined in table 17(b). The 

small value obtained for the coefficient of D 2 (0.47) (rather than 
the zero theoretical value) is due to imperfect fit. 

DISCUSSION 


Campari son of Fourier Transform Method, Derivative Method, 
and Prony*s Method of Obtaining a First Approximation 


The Fourier transform method appears to be very convenient, espe- 
cially if a harmonic analyzer is available. The derivative method does 
not appear to be practical unless the necessary higher derivatives are 
measured during the flight test. Although Prony*s method is limited to 
the case where the response is a sum of exponentials (the step response 
is included), it can always be used with a pulse input by fitting tbe 
free oscillations following the pulse. This determines the characteris- 
tic roots of the system, or the transfer coefficients in the den om i na tor 
of the transfer function. It is shown in reference 5 that once the 
transfer coefficients of the denominator have been determined, the trans- 
fer coefficients of the numerator of the transfer function can be com- 
puted by a linear calculation. In the case of the response due to a step 
input, the complete transfer function can be determined, as shown in the 
example, but this is a special case. 
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Determination of Lift Derivatives 


It can "be shown that, given a complete set of transfer coefficients 
(e.g. , those corresponding to q/8 and a/s) , the lift derivative 
is given by: 



b C, C 0 
iq u q 






m ¥ 0 



If Lg is negligible this reduces to 


iq_ 

mV 0 


= £2 a. 


The value of Iq, can also be determined from dynamic measurements 
of a single response at two center-of-gravity positions as follows: 


mV 0 Iy Iy 

Differentiating this with respect to x (assuming Iy and Mq 
do not vary appreciably with x) gives 


from which 


& = Ztt 

dx Iy 




The above methods of determining Iq, can be used no matter how 
large the value of Lg is. If Lg is small enough, the value of L a 
can be determined from static-stability measurements or from steady- 
turn measurements. (See subsequent section entitled "Relation Between 
Static and Dynamic Tests.") 


Basic Limitations in Determination of Moment Derivatives 


Various methods of computing transfer coefficients and lift deriva- 
tives are mentioned above; the computation of moment derivatives is not 
explained. For the two degree of freedom case discussed above, the 
unique determination of the damping moment derivatives from the airp lan e 
response to control deflections is impossible, as will now be shown. 
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If a ,, Da , q, 8, and Dq are known at four different instants 
of a transient disturbance, the moment equation can he used to form four 
simultaneous linear equations for the determination of the four unknowns, 
^Da* % %* Tlie coefficients of these unknowns are the 

values of a. Da, q and 8 at these instants. Unfortunately, these_ 
four values of a. Da, q and 8 are linearly dependent according to 
the lift equation, which states that 


aL^ + Da mY 0 + 8Lg — q mY 0 = 0 

Therefore, the four simultaneous equations formed to determine the 
moment derivatives are linearly dependent and, as shown in reference 7, 
cannot he solved uniquely for the unknowns. If the value of any one of 
the moment derivatives is known from other sources, then the values of 
the other three may he uniquely determined from the flight measurements. 
Otherwise, only combinations of the moment derivatives such as the 
following can he evaluated from the flight data: 

mV o + % 
mV 0 M 5 + Lg 

In practice, this indeterminacy prevents only the separation of the 
damping derivatives in the first of the above expressions, because they 
are of the same order of magnitude, at least for most configurations. 

The second of the above expressions often may he considered as 
practically equal to mY 0 for sufficiently high static margins, 
altitude, and wing loading. This can he more readily shown by taking 
the ratio of mYo^ to Mql^ which can he shown to he 

m ^oMq _ 2|iX 

\ L a Cm qc 
zV 

where 

H = -2L = 26.1 ifi.u 

ps„c V Po/ ° 


x = static margin 
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For w/Sy. = 100, c = 2, p/p Q = 0*53 (20,000 feet altitude), x = 0.2 
of M.A.C., and Cm^ c = —15 (conventional configuration), the ratio is 

w 

approximately 20. This means that the error in Mq, involved in 
neglecting the Mq_ term in the expression ^YqM^ + Mgl^ is 5 percent. 

The last of the above expressions, namely mYoMg + MqLg, also can 
he assumed equal to mVoMg to about 5 percent for conventional plan 
forms, so that Mg may be known to about this accuracy. 

If the tail length 2 defined as the distance from the airplane 
center of gravity to the center of pressure due to 5 is known, then 
Mg also can be found from Lg by the simple relation: 

Mg = — 2 Lg 

1 This additional relation between the derivatives Mg and Lg makes 
possible, in principle, the unique determination of all the derivatives; 
however, in practice, the uncertainty (especially for unconventional 
configurations) in the value of 2 and the low accuracy cf determining 
Lg (for conventional configurations where Lg is almost negligible) 
make the accuracy of the determination of the individual derivatives by 
this method very much lower than that for the transfer coefficients or 
combinations of derivatives. 

The indeterminacy of the moment derivatives from control— response 
data persists, in general, even when changes in forward speed are taken 
into account. Measurements of airplane response in the three longitu- 
dinal degrees of freedom (e.g., Y, n, and q) to a known control deflec- 
tion will yield values of 


Iry — M y 1^ 


^v Mdcg “ M v mV o 


and so forth. 

The derivative M y , a measure of the effect of airspeed changes 
alone on pitching moment, is caused primarily by power and Mach number 
effects. In cases where this derivative can be assumed negligible, 
measurements involving forward speed permit evaluation of all deriva- 
tives. It is interesting to note that a similar type of indeterminacy 
exists in the lateral motions if product-of— inertia terms are present. 
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Information Obtained From Tail Load Measurements 


If the loads on the tail of a conventional— - type airplane are 
measured during a maneuver, the value of downwash slope e a may be 
inferred, provided the following assumption is made: 

^ = aLb a t (1 " + 33(1 ^ % + 4 S 4. + SL 5 

By measuring the time histories of a (or n), 5, q and Lfc, and 
knowing 2 and V the values of , e a , and L§ nay be deter— 

“t 

mined using the same numerical methods outlined above. 

Use of If onaer odynaml c Forcing Functions 

It has been shown above that the indeterminacy that exists when the 
response to a control deflection is measured would be overcome if the 
location of the center of pressure of the lift due to control deflection 
were known. This immediately suggests a testing technique in which the 
force is applied at a known position in the aircraft by, for example, 
dropping a bomb or firing a gun or Jet normal to the longitudinal axis. 

If the known applied lift L acts at a known distance 2* behind 
the center of gravity, then the applied moment M is equal to — 2'L 
and the equations of motion are 


-ol^ —Da mY 0 + q mY 0 = L 


-oMq^ -Da -q Mg + Dq Iy = -2*L 


Solving these equations for a and q gives 


a 

L 


q. 

L 


1 D+ Mg — I’mYp 


MY, 




D 2 + b D + k 


+2 t mYo 1^3, + 2*1* 

XjfflVo 

D 2 + b D + k 


I y mY 0 
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where b and k are the same as before. The above transfer coeffi- 
cients (determined by the methods of the earlier part of this report) 
can obviously be broken down immediately to give and The 

determination of 1^ and Mq, from the transfer coefficients is also 
feasible, although somewhat less straightforward. 

If the bomb in the previous case is replaced by a gun firing 
vertically, the same method of analysis applies. If the impulse supplied 
by the gun is I (= mass of shell times muzzle velocity), the equations 
of motion are: 


_<xL — m\T Da + mY n q. = ID 
a o ° 


-aMa — Mj-^ Da — Mg q. + Iy. Dq = I 2'D 


where l * is now the distance between gun and center of gravity. It 
is not necessary to know the value of the impulse I in order to 
determine the derivatives — a fact which is of importance since the 
impulse is not as easy to' measure as the weight of the bomb in the 
previous example. 


Relation Between Static and Dynamic Tests 


In static or maneuvering stability tests, the rate of change of 
trim with control deflection is measured. This is, of course, the ratio 
of the constants in the transfer function. For example, the longitu- 
dinal maneuvering stability is determined by the normal acceleration per 
degree elevator deflection in a steady turn which is 

= ^ = Y 0 (^L S -Mglfr) 

D=0 k g^Y^ +1^1^) 

The variation of this with center-of-gravity position is 

d / n\ _ / n\ dk/dx 
dx \ ^ _ \ 8/ k 



n/s k 


The ratio 
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is the distance between the center of gravity and maneuver point (i.e., 
the point where k=0) . 

If the angle-of-attack change (as well as the normal acceleration 
change) in a steady turn is measured, we get 


/ a\ _ M5 mV 0 + Mq 

V 5 ^=0 ^ A k 


The variation with center-of-gravity position is 



The ratio 


a/S 

J“(a/6) 

dx 


does not accurately give the maneuver point unless L§ is negligible. 
If Lq is not negligible, as for a tailless airplane, it is necessary 
to measure 6 /a at a number of center-of-gravity positions to deter- 
mine the position for which 5 /a Is zero, in order to find the 
maneuver point. Thus, normal acceleration measurements enable a some- 
what more convenient maneuver point determination than angle-of-attack 
measurements. The ratio of n/6 to a/6 gives 

E = _ Z° ~ Iq, 

a g mV 0 M§ + Mq Lg 


If I 5 is negligible, this equals I^/mg. If L5 is not negligible, 
it is not possible to determine L a from steady-turn measurements. 

In static longitudinal-stability tests the values of V/S and Its 
variation with center-of-gravity position similarly determine the 
neutral— point position. It can be shown that 


V _ m 5 ^ ” % L S 

8 (ayf/V 0 ) - m v i<x. 
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and 

7/5 2% W - M T Ig 7 p 

jL(y/8) 2 W la 

CLX 

which is equal to the distance "between the neutral point and the center 
of gravity. The measurement of a/5 gives 


a _ mV 0 M§ + Mg L 5 
8 " -la Mg_ - mYoHx 


Just as in the case of steady turns, the value of 

g/5 

d/dx (a/5) 

does not quite give the neutral— point position unless Lg is negli- 
gible. Measurement of 7/8 and a /8 gives 

7 _ Iq M & - ^ Lg 
a- L 5 M v - Mg 2W/Y 0 

If Lg is negligible, the above expression reduces to 


2W 

but no static or steady— turn measurement at a single c ent er-of— gravity 
position can determine L or Lg if Lg is appreciable compared with 

V 


Aerodynamic lag 


The lag between aerodynamic forces and moments and angle of attack, 
control deflection, or pitching velocity is variously labeled as: 

Unsteady lift 
Indie ial lift 

Nonstationary or unstationary lift 
Oscillating airfoil effect 
Wagner effect 
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The assumption was made above that the aerodynamic forces are 
represented by 


Xi = Olfc + &Lg 


M = ai^ + Da + SMg + qM^ 


These equations indicate that the lift due to a sudden change in control 
deflection, angle of attack, or pitching deflection, and the pitching 
moment due to control deflection or pitching velocity develop instan- 
taneously. The lag in pitching moment due to angle of attack, supposedly 
represented by Md<x, is obviously an approximation to the relation 
between angle of attack and pitching moment due to the complex transient 
downwash disturbance at the tail. Actually the lift due to sudden 
changes in angle of attack, etc., lags somewhat behind the angle of 
attack at first. The same effect can be expressed in terms of the fre- 
quency response of lift due to sinusoidal variation in angle of attack 
and, therefore, is also referred to as the "oscillating airfoil" effect. 

The aerodynamic lag effect can be expressed in a number of ways, 
namely, 

(1) The expressions for L and M given above require addi- 
tional derivative terms on both sides of the equation to take account 
of the aerodynamic lag. The correct expression for L, for example, 
would be of the form: 

(a^ D n + a n _ i tf*- 1 + . . . a x D + 1)L = 

1^(1 + b 1 D + . . . bj^ D^)a + Lg(l + c-j_ D + » » . Cp D )5 

where m and r are less than n. 

(2) The relation between an aerodynamic force and airplane motion 
in a particular degree of freedom can be expressed as an "aerodynamic 
transfer function" which can be described in terms of an "indicial lift 
response" or an "oscillating lift response." The indicial lift response 
is the variation of lift with time for step function displacements in 
the partic ular degree of freedom. The oscillating lift response is the 
amplitude and phase of the lift with respect to displacement for 
sinusoidal motion in the particular degree of freedom as a function of 
frequency. 

The two methods of expressing aerodynamic lag are, of course, 
related by the fact that the expression given under (1) is the rational 
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fraction approximation to the "operational equivalent" of the indicial 
lift response mentioned under (2)« This means that the indicial lift 
response or oscillating lift response as calculated from (1) would he 
the same as that defined in (2). 

It is theoretically impossible (in the most general case) to deter- 
mine these unsteady lift functions from the response to control deflec- 
tion alone, as pointed out above for the determination of moment 
derivatives. Their importance in airplane maneuvers of the type dis- 
cussed in this report is believed to be small, however, except possibly 
in the transonic speed range. 


Suggestions for Future Work 


It would be of great interest to compute the error ratios of the 
parameters involved in the longitudinal and lateral dynamic airplane 
responses, for various forcing functions and airplane configurations. 
This would show what accuracy can be expected from various flight- 
testing techniques. The problem of taking account of errors in measure- 
ment of the forcing function (e.g., control angle deflection) remains 
to be solved. 

In applying the methods of this report it is necessary to assume 
the form of the transfer function in advance. Further research is 
desirable on a method of solution which would yield the form of the 
transfer function as well as the values of the transfer coefficients. 


COHCLUDUTG REMAKES 


A mathematical study of the determination of dynamic airplane 
parameters from flight measurements showed, as would be expected, that 
the determinability and accuracy of the parameters depend on many 
factors, such as type of input disturbance, the amount of response data, 
and the method of analysis. The determination of stability derivatives 
as distinct from transfer-function coefficients requires, in general, a 
more complicated flight— testing technique. In the case of longitudinal 
motions, for example, measurement of the response of the airplane to 
elevator motions theoretically will not enable derivatives to be 
uniquely determined; instead, the values of linear combinations of pairs 
are obtained. Practically, the static derivatives can be determined 
with reasonable accuracy in many cases. Measurement of the response to 
a known force, as applied for example by a bomb drop or gun recoil, 
yields additional data from which all the derivatives can be determined. 

Application of these analysis methods and principles should assist 
in developing testing and analysis techniques to satisfy the requirements 
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of various specific research programs. Most previous work concentrated 
on sinusoidal control-deflection inputs or converted data from the 
response to step or pulse inputs to frequency— response form for further 
analysis. Conversion of the transient response to the frequency 
response is not the only way to analyze the data to determine dynamic 
parameters. Methods such as Prony's and the "derivative method" are 
also available. Prony's method, in particular, appears to he a simple 
and direct method of determining the characteristic roots from free— 
oscillation data, and is also applicable to the response to a step input. 

In the determination of the dynamic parameters from flight data a 
rational least— squares procedure should be used to give the best fit to 
redundant data. Application of this procedure also gives an evaluation 
of the accuracy of a parameter, relative to that of the basic flight 
data. 


The methods of analysis presented in this report are applicable to 
more complicated problems than the ones used for illustration. The 
inclusion of additional degrees of freedom, and of additional higher- 
order derivatives to take account of unsteady lift effects is possible 
by the methods of this report, and remains a field for future work. 


Ames Aeronautical Laboratory, 

National Advisory Committee 
Moffett Field, Calif. , 


for Aeronautics, 
Dec. 19, 1950. 
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APPENDIX 


Details of Calculation of Parameters From First Approximation 


In the case where one transient response is measured (say q) the 
equations for the increments in the parameters become ; 


s iJ x j - s i 


1 = 1, 2, 3, k 


Sl3 = / d^)(sr) 4t Sl=_ / sf dt 


where 


^ .>»* r s( T,e-^ T aT + ° ^y < s .v r b(T ) 

Aj— X 2 J 0 X 2 —X 1 J 0 


C Ao+C r 


e 2 dT 


A x and X 2 are the roots of A 2 + boX + k Q = 0 


where 


are the first approximation 


= b— b Q 
X2 = k-k Q 

T <> = %- 0l 9o 

X * = S 0 


If X 1 and X 2 , the conjugate complex roots, are equal to 2 ± IV , then 1 


n =®I2 I (g z*cos Z*t+C, 2 sin 2't+C 0n sin Vt) [ &(T)e ^cos 2*TdT_ 

V l\ ^o ^o 1o / J 0 


x The symbols 2 and 2’ used in the appendix should not be confused with 
those in the text. 
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( 


‘4o 


C ln 2 cos 2*t + Co n cos 2't — 2’sin 2*t^ f 8(r)e ‘''sin 2‘Td.T 


— 2t 


/Jo 


■} 


In order to obtain the required derivatives and ^k it is 

n 9k 

OQ- dq r 

first necessary to obtain the derivatives —2 and — - and then use 


the expressions 




9X = 


&k = 

X j 

, X 2 

9b 


^Xg X 2"“Xg 

5k = 

9qc 1 

+ ^ 1 

9k 

^X^ Xg“"X^ 

^Xg X 2"“Xg 


The values of the partial derivatives are: 

t 


d^c = °1q X S + c o Q 

(X^^- -Ip) 2 


X_t 


/ &( T)e 


"*'2 T d T -e X i t 


J 8(r) e - A 'i T dT] 


Sqc 


9X, 


s expression obtained by interchanging X x and X 2 in above 


_ _X 

9c 


= e Xlt f 6(T)e _> ‘ lT 4T + I* 6 (r)e 

X-L-X 2 Jq Xg— A.J Jo 


-X-T , 

E dT 


9^ J* 8(T) e - X i T dT t e^ jf B(T) 

dC, 


— X P T 
e 2 dT 


■'OqL 


X — X 2 


Xp— X ■ 


Where 9 , q, and Dq are measured in a transcient response, the 
linearized equations for the increments become: 



J=1 


x 3 


- s i 


i = 1, 2, 3, ^ 
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«-» / ©© 



oT pin n T n T ^ 

51=5,9 ^ {0m ^ c) 3^ X («.-*.) dt+W B, X (E ^o) jjJ 


dt 


x a = Ab 


where 


0 C = — ^ 


x 2 = Ak 


x 3 = AC, 


x 4 - AC o q 


y" 1 s(T)dT 


c x q + Co X t nt -X T C x X + c o a 

-ter el x 8(T)e 4T + — — ^ 


^2(^2 ^*x) 


x 2 t r* 


J' S(T)e ^ 2>r dT 


C i a + c o, 


XXP 2 .Mr « T)t v tT + !h°X ,M r B(T) .-sv 

A 1 A 2 '“'O A 2 — ^1 '-'q 


In = 


DoLc = C, 8(t) + 


So*" 


i^o) 


x,t 


Ai 1 *~A<p 


X 


(c, X_ + Co. ]A 2 

r> / — \ — X,T V 2 ‘W 

6(T)e 1 dT + — - A 

V*i 


x„t 


/ 5(T)e X '2 t dT 


or if X x and X 2 are the conjugate complex roots equal to Z ± iZ* 
then 

c Z*k 


\ Sq 2 * cos in + c °io l sin 2,t + Cl q k sin v *) c 


-Cx k cos I«t - C 0 Z cos Z»t - C Q Z *sin Z*t + 

^0 'Iq ^o 



J °lo 


5(t) d T 
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.2!! iy c 

v [A 1 ? 


2* cos 2*t + C, 2 sin 2’t + Cr> sin 2*t } C - 


) 


(%c 


2 cos 2*t + Cr, cos 2*t — C-, 2* sin 2’t 


>] 


, 


It r r 

V { [ 2 °lc 


2 2*cos 2 f t + Cr, 2* cos 2 f t +[ C, 2 + C n „ 2 - 


(°* 


C, 2 * 2 ^ sin 2*tl C + rr 2 C. 2Z f + C 0rt v\ sin 2* t 

S 0 ) J LV *0 ) 

( C, 2* 2 - Co 2 - C, l 2 \ cos 2*tl dL + C n S(t) 

V V, ° 1 o ) J J V 


•where 


- / 8(t,s- It 


cos Z'TdT and D 


-f 


-2 T 

e sin 2' 


If the frequency response q = E sin (cot +9 ) to S = sin cot 

measured, the equations for the increments in h, k, Ci , and C 0 

* 

the first approximation are denoted "by: 

4 

i - 1, 2, 3, 4 


^ s iJ x j = s i 


J-l 


where 


•*<- -i(S © • - 1(5) © 


r— t ^Cp 

(^ - E c) ^ + “cp ^ (9m “ 9c) sf 
where the s umma tion is to he taken over the frequencies. 

1 a_ h„o) 


Coq 2 % 2 0)2 C 0 *o - “> 2 

up- Hq «,„;i 4-ot* m — *Q 


S c = 


— — = and tancp c = 

2.2 2 2 TC C, W 

(V® ) + “ 1+ 




T dT 


is 

over 
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TABLE I.- LEAST SQUARES CALCULATION OF FIRST APPROXIMATION TO TRANSFER 
COEFFICIENTS FROM FREQUENCY RESPONSE IN PITCHING 


a 

B 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

Row 

a> 

B 

<P 

(deg) 

sincp 

cosq> 

B sincp = 
B 

E coscp = 
A 

2 

A 

AB fu 

Ben 

2 2 
B o> 

2 

<0 

A 2 2 

A o> 

2 

Aon 

3 

Ben 

2 

B 

M 


8.8587 

184.8 

-0.06662 

-0.99778 

-0.5902 

- 8.8390 

78.1279 

5.2168 

-0.5902 

0.34834 

1 

78.1279 

-8.8390 

-0.5902 

0.34834 



10.0211 

183.4 

-.05878 

-.99827 

-.5890 

-10.0038 

100.0760 

11.7845 

-1.1780 

1.38768 

4 

400.304 

-40.0152 

-4.712 

.34692 

3 

3 

11.3560 

177.8 

.03847 

-.99926 

.4369 

-11.3476 

128.7680 

-14.8733 

1.3107 

1.71793 

9 

1158.912 

-102.1284 

11-7963 

.19088 

4 

4 

12.2791 

/ 

168.7 

.19543 

-.98072 

2.3997 

-12.0424 

145.0194 

-115.5926 

9.5988 

92.1370 

16 

2320.31 

-192.6784 

153.5808 

5.75856 

5 

5 

12.4966 

158.4 

.36748 

-.93003 

4.5923 

-11.6222 

135.0755 

-266.8631 

22.964 

527.230 

25 

3376.888 

-290.555 

574.0375 

21.08922 

6 

6 

12.0825 

148.8 

.51818 

-.85527 

6.2609 

-10.3338 

106.7874 

-388.1933 

32.5654 

1411.16 

36 

3844.346 

-372.0168 

1352.354 

39.19887 

7 

7 

11.3058 

140.6 

.63500 

-.77251 

7.1792 

- 8.7338 

76.2793 

-438.9119 

50.2544 

2525.50 

49 

3737.686 

-427.9562 

2462.47 

51.54091 

8 

8 

10.4232 

133.9 

.72007 

-.69390 

7.5054 

- 7.2327 

52.3119 

-434.2745 

60.0432 

3605.19 

64 

3347.962 

- 462.8928 

3842.76 

56.33103 

9 

9 

9.5589 

128.6 

.78119 

-.62429 

7.4673 

- 5.9675 

35.6m 

-401.0500 

67.2057 

4516.61 

a 

2884.499 

- 483.3675 

5443.662 

55.76057 

10 

10 

8.7671 

124.4 

.82551 

-.56439 

7.2373 

- 4.9481 

24.4837 

-358.1088 

72.3730 

5237.85 

100 

2448.37 

- 494.810 

7237.30 

52.3785 

11 


£ 





-91.0709 

882.5402 

-2400.8662 

314.5445 

17919.13 

385 

23597.40 

-2875.259 

21072.66 

282.9438 


Heal Equation o£ Condition 
Bm 1) — A4 + C 0 + Auf = Eg 


Imaginary Equation of Condition 
‘ ’ ~ % 


Ao> fc + Bk — -Ben 


Normal Equations Obtained by Minimizing E Eg 2 + E E-j- 2 


© + £©]b 

+ 

E ©c 0 - 

E ©Ci q = 

0 

[£ @ + £®]4 

- 

1 

V 

© 

w 

2 © c lq « 


£ © b - 2 ©4 

+ 

1°. C 0q 

= 

-£© 

E© 6 + E@4 


- 

S @Clq = 


41516.53 6 

+ 

314.5445 C 0q + 

2875.259 Ciq = 

0 

1165.484 4 

+ 

91.0709 C D - 

314.5445 Ciq = 

41516.53 

314.5445 6 + 91.0709 4 

+ 

10. c 

°q 

= 

2875.259 

- 2875.259 b + 314.5445 4 



385. C lq = 

21072.66 


b=8,0W5 


k=28.6l 6 


C 0q =-226.2 


C lg =^lM 


U) 
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TABLE 


II.- APPLICATION OP FOURIER TRANSFORM TO CONVERT THE RESPONSE TO 
A STEP INFOT INTO THE FREQUENCY RESPONSE 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

6* 

25 

26 

27 

20 

29 

30 

Row 

t 

q/6 

2t 

3t 

4t 

5t 

6t 

7t 

8t 

9t 

10t 

.la© 

•In @ 

•In ^6) 

ala© 

■10® 

■“© 

•In (S) 

■la © 

•In @ 

© 

m 

0 

eo# (?) 

Q 

s 

COB (l) 

COB (0) 

coa (?) 

cos (5) 

COB 

COB (l|) 

1 

0 

-0.00 

0 

0 

0 

a 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1,00000 

1.00000 

1.00000 

2 

.1 

"7.03 

.2 

•3 

.4 

.5 

.6 

.7 

.8 

•9 

1.0 

.19867 

.29552 

.38942 

.47943 

■36462 

.64422 

.71736 

.78333 

.64147 

.98007 

.95534 

.92106 

.87750 

.02534 

-76484 

.69671 

.62161 

.54030 

3 

.2 

-10.42 

.4 

.6 

.8 

1.0 

1.2 

1.4 

1.6 

1.8 

2.0 

.36942 

.56462 

.71736 

.04147 

.93204 

.98543 

.99957 

-97385 

.90930 

.92106 

.82534 

.69671 

.34030 

.36236 

.16997 

-.02920 

-.22720 

-.41615 

4 

-3 

-11.54 

.6 

•9 

1.2 

1.5 

1.8 

2.1 

2.4 

2.7 

3-0 

.56462 

.78333 

.93204 

-997*9 

■97385 

.66321 

.6754 6 

.*8738 

.14112 

.82334 

.62161 

■36236 

.07074 

-.22720 

-.50485 

-.73739 

-.90407 

-■96999 

5 

.4 

-11.48 

.8 

1.2 

1.6 

2.0 

2.4 

2.8 

3.2 

3.6 

4.0 

-71736 

-93204 

•99957 

■90930 

.67546 

.33*99 

-•03837 

-44233 

-.75680 

.69671 

.36836 

-.02920 

-.41613 

-.73739 

-.94222 

-99889 

-.09676 

-.63364 

6 

■5 

-10.00 

1.0 

1.5 

2.0 

2.5 

3.0 

3.5 

4.0 

4.5 

5.0 

.84147 

•99749 

.90930 

.59847 

.14112 

--33078 

-.73600 

-.97733 

-.958 98 

•54030 

.0707* 

-.41615 

-.80114 

-.98999 

-93646 

-.65364 

-.21080 

.28366 

7 

.6 

-10.24 

1.2 

1.8 

2.4 

3-0 

3.6 

4.2 

4.8 

5.4 

6.0 

.93204 

-97385 

.67546 

.14112 

-.44252 

-.87138 

-.99616 

-.77276 

-.27942 

■36836 

-.22720 

-.73739 

-.98999 

-.89676 

-.49026 

.08750 

.63469 

.96017 

8 

■7 

-9.66 

1.4 

2.1 

2.8 

3-5 

4.2 

4.9 

5.6 

6.3 

7.0 

.98343 

.86321 

.33699 

-.35078 

-.87158 

-.98243 

-.63127 

.01681 

-65699 

.169917 

-■50*85 

-.94222 

-.936*6 

-.49026 

.18651 

■77557 

.99986 

-73390 

9 

.8 

-9-18 

1.6 

2.4 

3.2 

4.0 

4.8 

5.6 

6.4 

7.2 

8.0 

-99957 

.67546 

-.05837 

-.73680 

-99616 

-.63127 

.11655 

•79367 

.98936 

-.02920 

-■73739 

-.99689 

-.6536* 

.08750 

•77557 

■993*8 

.60835 

-.14330 

W 

.9 

-8.90 

1.8 

.2.7 

3.6 

4.5 

5-4 

6.3 

7-2 

8.1 

9.0 

•973B5 

.*8738 

-.44252 

-97753 

-.77276 

.01681 

.7 9367 

■9698 9 

.*1818 

-.22720 

-.90407 

-.09676 

-.21080 

.63469 

•99986 

.60835 

-.8*33* 

-.91113 

11 

1.0 

-8.79 

2.0 

3-0 

4.0 

5-0 

6.0 

7-0 

8.0 

9-0 

10.0 

-90930 

-14112 

-.75680 

-.95898 

—879*2 

.65699 

.96936 

.41212 

—54403 

-.41615 

-.98999 

-.65364 

.28366 

.96017 

■75390' 

-.14550 

-.91113 

-.83907 

12 

1.1 

-6.70 

2.2 

3.3 

4.4 

5-5 

6.6 

7-7 

8.8 

9.9 

11.0 

.80650 

-•*577* 

-.951& 

—70554 

.31154 

.98817 

.58*98 

-.*575* 

-.99999 

-.58830 

-.90748 

-■30733 

.70067 

.95023 

•15337 

-.81109 

-.68919 

.00426 

13 

1.2 

-8.71 

2.4 

3-6 

4.8 

6,0 

7.2 

8.4 

9.6 

10.8 

12.0 

.67346 

-.44252 

-.99616 

-.27942 

-79367 

.85460 

-.17433 

—98093 

! —»53657 

-.73739 

-.09676 

.08750 

.96017 

.03035 

-.51929 

-.90469 

-.19433 

1 ,84305 

1* 

1-3 

-8.73 

2.6 

3-9 

5.8 

6.5 

7.8 

9.1 

ia.4 

11.7 

13.0 

-51550 

-.68777 

-.92581 

.21512 

•9985* 

.31910 

-.82783 

-.76198 

,42017 

-.85689 

-.72593 

.37798 

•97659 

.05396 

-.94772 

-.56098 

.64760 

•90743 

15 

1.4 

-8,77 

2.8 

4.2 

5.6 

7.0 

8.4 

9.8 

11.2 

12.6 

lfc.Q 

.33499 

-.87158 

-.63127 

.63699 

.85460 

-.36648 

-97918 

.03368 

.99061 

-.94222 

-.49026 

.77557 

•75390 

-.51929 

-.93043 

.20300 

.99943 

■13674 

16 

1-5 

-8.79 

3.0 

4.5 

6.0 

7.5 

9-0 

10.5 

12.0 

13.5 

1?.Q 

.14112 

-.97753 

-.279*8 

.93800 

.41212 

-.87970 

“.53637 

.80379 

.63089 

-■98999 

-.21080 

.96017 

•34664 

-.91113 

-.47554 

.84383 

.59492 

-73969 

17 

1.6 

-8.77 

3.2 

4.8 

6.4 

8.0 

9.6 

11.2 

12.8 

14,4 

16.0 

-.05837 

-.99616 

.11635 

.90936 

-•17433 

-.97918 

.23151 

.96566 

-.28790 

-.99829 

.08750 

.99318 

-.14550 

-.98*69 

.20300 

•97883 

-.25982 

-.93766 



* 


NACA TN 2340 




OS 


TABLE II.- CONCLUDED 


s 

> 

H3 

is! 

ro 

uo 


1 

2 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

Bow 


©*© 

©*© 

©*© 

©*© 

©«© 

®x@ 

©X® 

(3>@ 

©X® 

©=<§) 

©X® 

©X© 

©X@ 

®X® 

®x@ 


©*(§) 


Simpson's 

factor 

1 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-0.08000 

-0.08000 

-0.08000 

-0.08000 

-0.08000 

-0.08000 

-0.08000 

-0.08000 

-0.08000 

1 

2 


.1 

-1.40062 

-2.0Q342 

-2.74541 

-3.37998 

-3.98057 

-4.54175 

-5.05739 

-5.52248 

-5-93236 

-6.90949 

-6.73515 

-6.*93*7 

-6.18694 

-5.81865 

-5.39212 

-4.91181 

-4.38235 

-3.80912 

4 

3 


.2 

-4.05776 

-5-88334 

-7.47489 

-8.76816 

-5.71186 

-10.26839 

-10.41552 

-10.14752 

-9.47491 

-9.59745 

-8.60004 

-7.25972 

-5-62993 

-3-77579 

-1.77109 

.30426 

2.36742 

4.33628 

2 

it 


.3 

-6.51571 

-9.03963 

-10.75574 

-u. 51103 

-II.23823 

-9.96144 

-7.79481 

-4.93197 

-1.62852 

-9.52442 

-r. 17338 

-4.18163 

-.81634 

2.62189 

5.82597 

8.50948 

10.43297 

11.42448 

4 

5 


.4 

-8.23529 

-10.69982 

-11.47506 

-10.43876 

-7.75428 

-3.84569 

.67009 

5.00013 

8.68806 

-7.99823 

-*.15989 

• 33522 

4.77740 

8.46524 

IO.81669 

11.46037 

10.29480 

7.50379 

2 

6 


.5 

-9.15519 

-10.85269 

-9.89318 

-6.51135 

-1.53539 

3*81649 

8.23390 

10.63553 

10.43305 

-5.67846 

-.76965 

4.52771 

8.71640 

10.77109 

10.18868 

7.IU60 

2.29350 

-3.00622 

4 

7 


.6 

-9.54409 

-9.97222 

-6.91671 

-I.44507 

4.53140 

8.92498 

10.20068 

7.91306 

2.86126 

-3.71057 

2.32653 

7.550 SI 

10.13745 

9-18282 

5.02026 

-.896OO 

-6**9923 

-9.83214 

2 

0 


.7 

-9.51945 

-8.33861 

-3.23600 

3.38853 

8.41946 

9.49047 

6.098O7 

-.16238 

-6.34652 

-1.64191 

4.87685 

9.10184 

9.04620 

4.73591 

-I.80169 

-7.49201 

-9.65865 

-7.28267 

4 

9 


.0 

-9.1T605 

-6.20072 

.5358* 

6.94742 

9.14475 

5.79506 

-1.06993 

-7.28589 

-9.08232 

.26806 

6.76924 

9.16430 

6.00042 

-.80325 

-7-11973 

-9.11739 

-5.58465 

1.33569 

2 

10 


.9 

-8.66727 

-3.80368 

3-930*3 

8.70002 

6.87756 

-.14961 

-7.06366 

-8.63202 

-3-66787 

2.02206 

8.04622 

7.98116 

1.07612 

-5.64874 

-8.89675 

-5-41432 

2.16751 

8.10906 

4 

n 


1.0 

-7.99275 

-1.24044 

6.65227 

8.42891 

2.45610 

-5.77494 

-8.69647 

-3.62253 

4.78194 

3.65796 

8.70201 

5.74550 

-2.49337 

-8. *3969 

-6.62678 

1-27895 

8.00883 

7.375*3 

2 

12 


1.1 

-7.03395 

1.37233 

8.27892 

6.13820 

-2.71040 

-8.59708 

-5.08880 

3.98060 

8.69991 

5.11995 

8.59108 

2.67377 

-6.165*3 

-8.26700 

-1.33*32 

7.05648 

7.73595 

-.03706 

4 

13 


1.2 

-5.88326 

3.85435 

8.67655 

2.43375 

-6.91207 

-7.44357 

1.51841 

0.54390 

4.67352 

6.42267 

7.81078 

-.76213 

-8:36308 

-5.29073 

4.52302 

6-77665 

I.6926I 

-7.3*993 

2 

14 


1.3 

-4.50032 

6.00423 

8.08232 

-I.87800 

-8.71725 

-2.78574 

7.22696 

6.65209 

-3.66808 

7.48065 

6.33737 

-3-29977 

-8.52563 

-.*7107 

8.27360 

*.69736 

-5.65355 

-7.92204 

4 

15 


1,4 

-2.93786 

7.6W 

5.53624 

-5-76180 

-7.49484 

3.21403 

8.58741 

-.29485 

-8.68765 

0.26327 

4.29958 

-6.80175 

-6.61170 

4.55417 

8.15987 

-1.70O31 

-8.76500 

-1.19921 

2 

16 


1.5 

-1.24044 

8.592*9 

2.45610 

-8.24502 

-3.62253 

7.13256 

4.71645 

-7 .06531 

-5.71605 

8.70201 

1.85293 

-8.43989 

-3*04697 

8.00883 

4.18000 

-7.41744 

-5.22935 

6.67768 

4 

17 


1.6 

.5U90 

8.73632 

-1.02214 

-8.67669 

1.52887 

8.58741 

-2.03034 

-0.46884 

-2.52488 

8.75500 

-.7673' 

-8.71019 

1.27604 

8.63573 

-1.78031 

-8.53172 

2.27862 

8.39860 

1 

18 

Simpson Summation 

-287. 27402 

-108.85646 

-25.45190 

-79.07855 

-95-98573 

-32.59403 

4.64220 

-28.28000 

-41.26100 

.76806 

93-55412 

14.63327 

-23.57994 

40.05311 

68.30965 

20.39870 

-3.94770 

28.95594 


19 

^ x row 18 

-9.57580 

-3.62055 

-.84840 

-2.63595 

-3.19952 

-1.08647 

.15*7* 

-.94267 

-1.37537 

.02560 

3.11847 

.48778 

-.78600 

1.33510 

2.27699 

".6ra?6 

-.13159 

.96520 


20 


■» 

2 

3 

4 

5 

6 

7 

8 

9 

10 

2 

3 

4 

5 

6 

7 

8 

9 

10 


21 

u x row 19 

-19.15160 

-10.08562 

-3.39360 

-13-17975 

-19-19712 

-7.60529 

1.23792 

-8.48403 

-13.75370 

.05120 

9.35541 

1.95112 

-3-93000 

8.01060 

15.93893 

5.43968 

-1.18431 

9.65200 


22 

cosine 1.6 m 

-.99829 

.08750 

•99318 

-14550 

-.98469 

.20300 

.97824 

-.25978 

-.95767 

-.05837 

-.99616 

.11655 

•96936 

-.17433 

—97918 

.23151 

.96566 

-.28790 

sin 1.6u 

23 

8.77 x row 2? 

-6.75500 

•76738 

8.71019 

-1.27604 

-8.63573 

1.78031 

8.57916 

-2.27827 

-8.'390T7 

.51190 

-8.73632 

1.02214 

8.6766 9 

-I.52887 

-8.58741 

2.03034 

8.46884 

-2.52488 

8.77 x row 22 







In Phase component 








uut or Phase Component, B 






row 21-row 23 

-10.3966 

-11.6530 

-12.10# 

-11.9037 

-10.5614 

-9*3856 

-7.3412 

-6.2058 

-5.3549 

-.4607 

.6190 

2.9733 

*.7*67 

6.4817 

7.3515 

7-4700 

7.2845 

7.1271 



Simplon Sumntlon «* y (no. frcra colu m ) X (corresponding Simpson's factor trcm column 49) 

tSo 


By Simpson's r^ilotho Integral "of colum n la given by ~ times the Slnpion'e Sunmtlon. 



UO 

VO 






























































TABLE III.- LEAST-SQUARES DETERMINATION OE FIRST APPROXIMATION TO TRANSFER 
COEFFICIENTS FROM TEE RESPONSE TO A UNIT STEP INPUT USING DERIVATIVE METHOD 



MCA : 
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TABLE IV.- APPLICATION OP FRONT'S METHOD TO HIE RESPONSE TO A STEP INPUT. 

(a) Computation of b and k. 




-7.05 -10.42 
-10.42 -11.54 
-11.54 -11.48 
-11.48 -10.88 
-10.88 -10.24 





0.0064 

49.7025 


-80.53 -89.35 


-9.18 104.8576 98.9184 94.0032 
-8.90 93.3156 88.6188 85.9740 
-8.79 84.2724 81.7020 80.&22 


-91.09 824.0673 832.3838 810.6928 


8 

9 

10 

© 

X 

© 

@* 

© x ® 

0.6336 

49.7025 

73.4610 

81.3570 

108,5764 

120.2468 

119.6216 

133.1716 

132.4792 

125.5552 

131.7904 

124.9024 

117-5552 

118.3744 

1U.4112 

105.1008 

104.8576 

98.9184 

94.0032 

93.3156 

88.6788 

85.9740 

84.2724 

81.7020 

80.6922 

79.2100 

78.2310 

810.6928 

903.2709 

910.0308 


Equation of Condition: 81^ + q m+1 + «a + " E 

Normal Equations 2 @ax+£©Ua+ 2(|)a a =-S(D 824.0673 ax + 832.3638 aa - 80.53 « 3 - -810.6928 

Obtained by Min- 2 © n x + £ ® a a + 2 ® a 3 * -S © 832.3638 a x + 903.2709 da - 89.35 83 = -910.0308 

imizing £B a 

2* © ax + 2 © da + E @ a 3 * — £© -80.53 ax — 89.35 82 + 9 as ■ 91.09 

a x = 0.432127 da =-1.227683 ag - 1.799523 


let x = e " then x 8 + dax + ax = 0 

Solution for x - xx ± ixa = 0.61382 ± 0.2352l4i = Be* 10 = 0.657364 e ±0 ’ 30s5al 


b = —(lx + X a )= — 


log Be 1 + . 


log B g loj 

0.1 C 


_ax . (0.4^127). 8.3904 


* „ Xxla = (l°g^)(lo g Be-4g) _ (log R) a ♦ f _ (log 0.6573$* ♦ (0.36558)* _ 

0.1 0.01 0.01 

Equation for 1 = 1 +il ' X 2 + bX + 4 = 0 

X = -4.1949 ± 3. 6618 i then I = -4.1949 
I* = 3.6618 


In addition 


q. 2 = -8.802034 

ax + a 3 + 1 






























I 


TABLE IV.- CONCLUDED 
(to) Computation of C lq and Cq^. 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 


t 

q. - a«o = 
q _ 8.80 

2 t 

e-l* 

© 

X 

© 

■9 

sin © 

cos 0 


©* 

© X ® 

© 

X 

© 

© X © 

■a 


8.72 

0 

1.00000 

8.72 

0 

0 

1.00000 

0 

1.00000 

0 

0 

8.72 

'3 


1.75 

-.41949 

1.52120 

2.66210 

.36618 

.35804 

.93371 

.128193 

.871814 

.334306 

.953138 

2.485629 

B 

(SR 

- 1.62 

-.83898 

2.31405 

- 3.74876 

.73236 

.66861 

.74361 

.447039 

.552956 

.497185 

-2.506458 

-2.787615 

il 


-2.74 

-1.25847 

3.52014 

- 9.64518 

1.09854 

.89053 

.45492 

.793044 

.206952 

.405120 

-8.589322 

-4.387785 

5 

9 

- 2.68 

-1.67796 

5-35484 

- 14.35097 

1.46472 

-99437 

.10592 

.988772 

.011219 

.105324 

-14.270174 

-1.520055 

6 

.5 

- 2.08 

-2.09745 

8.14578 

- 16.94322 

1.83090 

.96638 

-.25713 

.933890 

.066116 

-.248485 

-16.373589 

4.356610 

7 

.6 

- 1.44 

-2.51694 

12-39137 

- 17-84357 

2.19708 

.81025 

-.58609 

.656505 

.343501 

-.474879 

- 14.457753 

10.457938 

8 

.7 

-.86 

-2.93643 

18.84033 

- 16.20268 

2.56326 

.54669 

-.83734 

.298870 

.701138 

-.457765 

-8.857843 

13.567152 


.8 

-.38 

-3.35592 

28.64560 

-10.88533 

2.92944 

.21064 

-.97756 

.044369 

.955624 

-.205913 

-2.292886 

10.641063 



-.10 

-3.77541 

43.59751 

- 4-39575 

3-29562 

-.15333 

-.98817 

.023510 

■976480 

.151516 

.674000 

4.343748 


9 

.01 

-4.19490 

66.35373 

.66354 

3.66180 

-.49705 

-.86772 

.247059 

.752938 

.431300 

-.329813 

-.575767 


£ 








4.561251 

6.438738 

.537709 

-66.050700 

45.300918 


Equation, of Condition q — q«> — e LZ (H sin l*t + N cos l *t) * E 


Normal Equations 
Obtained by Minimizing EE 2 


S © M + £ ©W =■ E © 


E@M + Z@R - £ © 

Response Equation q = - 8.802 + er “ 8B 6 (- 15.46 sin 3 . 662 t + 8.327 cos 3 - 662 t) 

5 662 

Laplace Transform of e-^* 195t sin 3 . 662 t - p 5 + bp t k 

Laplace Transform of e ~ 4 - l95 t cos 3 - 662 t = pi + ' bp^^k 

Transfer Function of i = - 8.802 - 15.46 - - + 8.327 

6 D®+ 8 . 39 D + 30.96 D 2 + 8 . 39 D + 30.96 

- 0 . 47 D 2 - 91 . 54 D - 258.4 
D e + 8 . 39 D + 30.96 


4-561251 M + 0 . 537709 H - - 66.0507 
0.537709 M + 6438738 N > 45.300918 
M - - 15.4625 N . 8.3270 


C lq = - 91.54 = - 258.4 


-t=" 

LO 


MCA TN 23^ 
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Figure /.- Frequency response of pitching velocity of example 

oscillations. 
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Figure 2 - Illustration of transient response to elevator pulse . 
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figure 5. 
ve/oci, 
input. 



